Eu^+ spin dynamics in the filled skutterudites EuM4Sbi2 (M=Fe, Ru, Os) 
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We report evidence for a close relation between the thermal activation of the rattling motion of 
the filler guest atoms, and inhomogeneous spin dynamics of the Eu^''" spins. The spin dynamics 
is probed directly by means of Eu^"*" electron spin resonance (ESR), performed in both X-band 
(« 9.4 GHz) and Q-band {^ 34 GHz) frequencies in the temperature interval 4.2 < T < 300 K. A 
comparative study with ESR measurements on the /J-EusGaieGeao clathrate compound is presented. 
Our results point to a correlation between the rattling motion and the spin dynamics which may be 
relevant for the general understanding of the dynamics of cage systems. 



I. INTRODUCTION 

The discovery of localized plionon modes in metallic 
compounds opens an avenue for the study of electron- 
phonon and phonon-phonon interactions in solids. The 
family of the filled skutterudite compounds are among 
the materials where this scenario is believed to take 
place^. These are cage systems inside which a guest, or 
filler, atom may perform relative large excursions. These 
excursions are described in terms of localized and isolated 
phonon modes (Einstein modes), usually called rattling 
modes. In this sense, the guest is referred as the rattler 
atom and its dynamics are fully characterized by a single 
parameter Oe, the Einstein temperature^. 

Filled skutterudite compounds have the general for- 
mula RT4X12, where usually R is a rare earth or acni- 
tide; T is a transition metal (Fe, Ru, Os or Co), and X 
is a pnictogen (P, As or Sb). These compounds crystal- 
lize in the LaFe4Pi2 structure with space group Imi and 
local point symmetry T/j for the R ions^. This structure 
hosts a wide range of physical properties, including exotic 
strongly correlated ground states'' as well as presenting 
a promising potential for application in the construction 
of thermoelectric devices^!. 

The latter feature is argued to be closely related to 
the dynamics of the rattler R atom^ . The rattling modes 
seem to promote the dampening of the thermal conduc- 
tivity mainly through the incoherent scattering of the 
Debye phonons. This leads to the concept of a phonon- 
glass type of heat conduction, which gives a fair picture 
of experimental results. Two important approximations 
to the rattler ion excursions are commonly adopted in 
this scenario: a description in terms of isolated (non cor- 
related) and localized (non dispersive) phonon modes^. 

Extensions of this scenario have been a matter of in- 
tense discussion in the field. One may find strong ex- 
perimental evidence for the dispersion of these modes'" * 
which raises the question of the validity of the indepen- 
dent rattler approximation. It also raises the question of 
weather or not an "electron glass", meaning a coupling 
between the guest ion dynamics and electronic degrees 
of freedom, could also be realized in these materials. In- 



deed, there are some recent theoretical proposals in this 
directiorPEl. In experiments, there is evidence for an 
interplay between the rattling motion and quadrupolar 
fluctuations with the unconventional superconductivity 
in PrOs4SbiJi2. Moreover, ^^^La nuclear magnetic reso- 
nance (NMR) show a correlation between the d band of 
transition metal and the rare earth rattling motion real- 
ized through electron-phonon coupling^'^. In this direc- 
tion, the rattling modes may also have an important role 
in the strongly correlated electronic phenomena found in 
cage systems. For applications, a better understanding 
of the description of the rattling modes may possibly lead 
to a more efficient design of thermoelectric materials. 

In a previous work, we presented electron spin reso- 
nance (ESR) as a probe for the guest ion dynamicJi^l 
We showed that the Yb'^^ spectra in Cei_a;Yb2;Fe4Pi2 
{x ^ 0.002) clearly reveal the R ion dynamical behav- 
ior. Taking advantage of this novel application of the 
ESR technique, here, we investigate the EuT4Sbi2 (T = 
Fe, Ru, Os) skutterudites, to observe whether the Eu^+ 
ESR spectra is a good probe for their rattling behavior. 
For comparison, we also present measurements on the (3- 
EugGaigGeso clathrate compound, another well-known 
cage system inside which Eu^+ behaves as a rattler ion. 

EuT4Sbi2 (T — Fe, Ru, Os) are metallic systems that 
undergo a ferromagnetic transition around Tc = 90 K, 
= 4 K and = 9 K, respectively !^. Close investigation of 
the T = Fe system indicates that, in fact, a ferrimag- 
netic transition takes place at Tc = 90 K, as a conse- 
quence of an antiferromagnetic coupling of the Eu and 
Fe moment&i^-^. The rattling behavior of the Eu^"*" ions 
was studied by structure refinement of x-ray diffraction 
measurements!^^ ^^-^^ g^Y^Q j^y extended x-ray fine struc- 
ture (EXAFS) measurements^^. The results obtained 
from both methods are in close agreement. The first 
measurements gave, for T = Fe, Ru, Os, 0e = 84, 78, 
74 K respectively, whereas EXAFS gave, for T=Ru, Os 
9e ~ 81, 78 K respectively. In both measurements, 
no signs of static off-center displacement of the Eu^^ 
were found. X-ray absorption near edge spectroscopy 
(XANES) and also susceptibility measurements, point to 
a predominant 2-\- valence state of the Eu atoms, al- 



though a slightly mixed valence state cannot be com- 
pletely excluded {v < 2.1)1^. In particular, for EuT4Sbi2 
(T=Ru, Os), the susceptibility measurements are com- 
patible with the full Eu^^moment in the temperature 
interval 2 < T < 400K. 

The type-I clathrate compound /^-EugGaigGeso is a 
metallic compound which undergoes a ferromagnetic 
transition at Tc — 35 K. In type-I clathrates, one finds 
two distinct types of cages that eventually lead to differ- 
ent dynamical behaviors of the guest elements. In partic- 
ular, in /S-EugGaigGeao, the majority of the guest Eu~+ 
ions are subjected to a potential with four symmetric 
off-center energy minima. In these cages, it was found 
that the dynamical behavior of the Eu^+ ions includes 
not only an off-center thermal activated rattling mode, 
but also quantum tunnel ing b etween the low temperature 
potential energy minimeP^ESl 

In this work, we report a non-trivial evolution of the 
X-band ESR spectra of EuT4Sbi2 (T = Ru, Os). We 
shall discuss that, in the low temperature region (4.2 < 
T < 150 K), the rattling motion couples with the Eu^+ 
spins. This gives rise to inhomogeneous spin dynamics, 
which manifests in the ESR linewidth (AH) behavior. 
We shall also discuss how our Q-band ESR spectra give 
support to these findings. 



II. EXPERIMENT 

For this work, we use single crystals of EuT4Sbi2 
(T=Fe, Ru, Os) grown in Sb-fiux as described in Ref.^'^. 
The resulting filled skutterudite structure was checked 
by x-ray powder diffraction. High filling rates were con- 
firmed by refinement methods. The ESR experiments 
were carried on crushed small pieces of single crystals 
with high grain size homogeneity. The ESR spectra were 
taken in a X-band and Q-band Bruker spectrometer us- 
ing appropriate resonators coupled to a T-controUer of 
helium gas flux system. The experiments covered the 
temperature interval 4.2 < T < 300 K. For compari- 
son, X-band ESR experiments were performed on a single 
crystal of the clathrate compound /S-EugGaigGeso. 

ESR detects the power P absorbed from the transverse 
magnetic microwave field as a function of the static mag- 
netic field H. The sensitive of the instrument is improved 
by applying a lock-in technique with field modulation. As 
a result, it is the absorption derivative jj^ which is ob- 
served . In our experiments, the ESR spectra showed a 
single Dysonian lineshape, described by equation: 
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this lineshape contains an a = D /A parameter express- 
ing the ratio between the dispersion {D) and absorption 
{A) of the microwave radiation when it probes a metal- 
lic materiaPD, This D /A parameter appears due to skin 
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Figure 1: X- band ESR spectra of EuT4Sbi2 (T = Ru, Os). 
Figures [if a)-(b) give an overview of the spectra evolution in 
the temperature interval 4.2 < T < 300 K and show in detail 
the spectrum for T = 4.2 K. As the temperature increases, 
the spectra broaden (Region 1) and then undergo a narrowing 
process (Region 11) , reaching temperatures (Region 111) where 
a Korringa Hke relaxation process describes the behavior of 
AH as T increases further. We also indicate the magnifying 
factor for each region, in the lower left corner of the figures. 



depth effects. A small excess of EuRu4Sbi2 was synthe- 
sized and the resulting crystals were crushed and sieved 
into fine powder. The powdered crystals were investi- 
gated in X-band and we obtained very similar results, 
but in this case the lineshape was completely symmetric 
{D = 0, no skin depth effects). This investigation was 
important rule out any surface effect in our results. 



III. RESULTS AND DISCUSSION 

Figures [l] [2] and [3] present the relevant experimental 
results for the skutterudite compounds and figure |4] gives 
a brief account of our results for the clathrate compound. 
In figures llTa) and flTb) , we give an overview of the X- 
band ESR spectra obtained for EuT4Sbi2 (T ~ Ru, Os). 
It is well seen that the thermal activation induces a fast 
broadening of the spectra until T sa 75 K, when a nar- 
rowing process sets in up to T « 150 K. Abo ve thi s T, 
the spectra broaden as in a Korringa relaxatiorpiEU 

As pointed out in these figures, there are 3 regions 
inside which the spectra evolve with distinct characteris- 
tics. The behavior in region III (T > 150 K) is expected 



to occur in any metal, but the behavior found in regions 
I and II is unique. For T = Fe, only the ferromagnetic 
modes (for T < T^, = 90 K) were observed for this com- 
pound (not shown). The non-observation of the ESR 
signal in EuFe4Sbi2 may be related to one or both of the 
following reasons: i) as presented in figures [21(a) and[2[b) 
(and also in^fa) and^Tb)) AiJ found in these compounds 
roughly scale with Tc {AHr^/AHqs ~ T^^'/Tg'). In 
a extrapolation, this would mean that AHpe would be 
of the order of several kOe, that would prevent its ob- 
servation; ii) EuFe4Sbi2 is a ferrimagnetic compouncjISl 
where the ordered Fe 3d orbitals partially compensates 
the Eu^+ moments. The effect of this local antiferromag- 
netic coupling between the Fe 3d moments and the Eu^"*" 
4/ moment may be a huge shift of the resonance that 
would also preclude its observation. 

Figures [2Fa) and |2jb) give a better view of the above 
cited regions by presenting the X-band ESR AH for both 
compounds. The behavior of AH resembles a coher- 
ence peak, as found in some NMR experiments22, when 
a phase transition is approached. From region I, for both 
compounds, as T increases, AH broadens by « 800 Oe. 
Following this broadening, around T « 75 K, both spec- 
tra gradually narrow by « 300 Oe (region II). At T « 150 
K, they broaden in a Korringa like relaxation as evi- 
denced by the fitting curve (region III). The measured 
Korringa rates are b = 5.5(1) Oe/K for T = Ru, and 
b = 6.3(3) Oe/K for T = Os. 

The insets in these figures show the normalized ESR 
intensity {xesr{T)/x{T = 300 K) for both compounds 
. A comparison with the Curie law shows that in all 
regions the spectra behave as a resonance arising from 
localized electronic states. This comparison also sug- 
gests that a small enhancement of xesr{T) is in or- 
der at low temperatures. This could be ascribed to the 
proximity of ferromagnetic transition. The ESR inten- 
sity was also normalized taking into account the varia- 
tion of the skin depth resulting from the T-dependency of 
resistivity^'^. In the experiment with the powdered sam- 
ple (EuRu4Sbi2), where the skin depth is bigger than 
the size of grains, xesr{T) also follows a Curie-like be- 
haviour. 

The similarities between the broadening and narrow- 
ing phenomena in regions I and II in both compounds 
indicate that these phenomena should have a common 
origin. In these compounds, the filler Eu atom is known 
to rattle with Einstein te mperat ures 0^; = 78 K (T = Ru) 
and 6'_E = 74 K (T = Os)Q^ESl. These temperatures are 
very similar to those where AH approaches its maximum 
before starting to narrow. We suggest that the thermal 
activation of the rattling motion gives rise to the behavior 
of AH in regions I and II. 

The thermal activation has its origin in the energy dis- 
persion of an Einstein (harmonic) oscillator E — {n + 
^)Huje {^e = kB^E-^ which, together with the proper 
statistical distribution, gives the number of active oscil- 
lators at a given temperature. Hence, in region I, there 
are rattling and non-rattling resonating ions. At low-T, 



2.0 



1.6 



1.2 



<] 0.8 
0.4-1 



T=Ru 




100 200 



0\ ^'^' 

■ Korringa law 

, (b = 5.5(1) Oe/K) 




(a) 




50 100 150 200 250 300 

T(K) 

Figure 2: X- band ESR linewidth of EuT4Sbi2 (T = Ru, 
Os). The same regions I, II and III discussed in fig[Ifa)-(b) 
are indicated in the above figures. These data show clearly 
the existence of such distinct regions. In the insets of figures 
[2fa)-(b), we present the X- band ESR intensity {xesr )■ It 
is shown that in all regions, xesr closely follows a Curie 
like behavior, which indicates that our resonance is due to 
localized electronic states. 



the lines are relatively narrower due to the spin-spin Ex- 
change interactioEp^. As T increases, the inhomogeneity 
of the spin dynamics, implied by the presence of rattling 
and non-rattling resonating centers, leads to the broaden- 
ing of AH. In region II, most of the ions are rattling, and 
an homogenization process ensues, leading to narrowing 
of AH. In region III, virtually all ions are rattling and 
the system has become quite homogeneous. 

As a plausible scenario for the origin of this correla- 
tion between the rattling modes and AH, we referred 
to the Ruderman-Kittel-Kasuya-Yosida (RKKY) inter- 
action, which gives the exchange coupling between the 
localized spins in rare earth metals^^. In skutterudites, 
the localized spins oscillate in their Einstein frequencies 
and this oscillation may give rise to random anisotropics 
in the exchange coupling J of the RKKY interaction, re- 
sulting in an inhomogeneous spin dynamics of the Eu^^ 
spins. Hence, in region I, the anisotropics of J broaden 
the resonance lines and prevent the narrowing of the spec- 
tra to occur. Once the rattling is fully activated, a nar- 
rowing process takes place, giving rise to the behavior of 
AH in region II. In region III, the narrowed spectra relax 
as in a Korringa- like relaxation. 

The above discussed behavior of the system is very sim- 
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Figure 3: Q- band ESR linewidth of EuT4Sbi2 (T = Ru, Os). 
Here, no clear relation with the rattling modes is observed. 
For T> 30 K, the spectra evolve as in a Korringa-like relax- 
ation and the slowing down of the relaxation is seem below 
this temperature. The onset of the later phenomenon starts 
at a slightly higher temperature in the case T=Os. . In 
the insets of figures [3ra)-(b), we present the Q-band ESR g- 
values. It is shown that at T ~ 35 K and T ~ 50, for T =Ru 
and T=Os, respectively, there is an increase of the g-values 
refiecting the development of a ferromagnetic internal field. 
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Figure 4: X- band ESR measurements of /3-EugGai6Ge3o. 
Figure p[a) gives an overview of the spectra in the investi- 
gated temperature interval. At high temperatures, the spec- 
tra is isotropic and the ESR intensity (not shown) follows a 
Curie law. As the magnetic transition is approached, strong 
anisotropy develops, as demonstrated by the two insets. Fig- 
ure p] (b), show that above T = 50 K, a Korringa like re- 
laxation dominates the behavior of Ai/, whereas at low T, 
A_H" broadens due to magnetic fluctuations. The inset shows 
strong anisotropy reflected in the g-values. 



ilar to a spin-glass, as probed by an ESR experiment^'* 
with a glass transition temperature Tq « 9e- For tem- 
peratures well below 6e, the rattling is not activated, and 
the low-T AiJ is only due to dipolar interactions and the 
spin-spin exchange narrowing. However, with increas- 
ing temperature, the spin-glass behavior manifests in the 
evolution of AH . In this sense, the peaks of AiJ, as 
shown in Fig. [2] are indeed coherence peaks, although 
we cannot claim that they express a phase transition to 
a true spin-glass state. 

In contrast to the above described unique behavior of 
the X-band spectra, the Q-band measurements, figures 
^a) and p[b) , do not present this conspicuous influence 
of the Eu^ rattling behavior. For T > 45 K, the spectra 
have a A_ff sa 250 Oe and evolves as in a Korringa-like 
relaxation with b = 6.7(1) Oe/K and h = 10.6(1) Oe/K 
for T — Ru and T — Os, respectively. In this T- region, 
the relatively narrow Aif allows for a precise determi- 
nation of the g-values. The results are g — 1.961(3), 
for T = Ru, and g = 1.92(1), for T =0s. In addi- 
tion, at low-T, both values increase, reflecting the in- 
ternal flelds of the low-T ferromagnetic state previously 



reported for the compounds^. These g-values reflect a 

huge g-shift (A5 = Qinsulator - 9experiment), which is duC 

to the Exchange interaction with conduction electrons (a 
polarization effect analogous to the Knight shift in NMR 
experiments), when compared to g = 1.993(2) found in 
insulatorgSH. This negative Ag is direct evidence for a co- 
valent mixing between the Eu^+ / orbitals and the con- 
duction band d orbitals. 

In a single band picture^, these Ag imply a Korringa 
rate of 6 « 24 Oe/K and 6 sa 124 Oe/K for T = Ru 
and T = Os, respectively. This strongly contrasts with 
the experimental results, reflecting the existence of a q- 
dependency of the exchange coupling J{q)'^^ . This is evi- 
dence that the Eu-Eu interaction is not due to a constant 
J exchange coupling, giving support to our claim of an 
inhomogeneous RKKY interaction. 

An alternative scenario for a correlation of the rattling 
behavior with the ESR spectra would be related to crys- 
tal field inhomogeneities. In this picture, the Eu^"*" while 
performing its excursion within the cage would experi- 
ence slightly different crystal field parameters and/or a 
crystal field environment of lower symmetry. This was 



shown to be the case in our previous work on the Yb'^"'" 
resonance^. In this work, however, the broadening of 
AiJ in regions I and II, instead of being enhanced, is 
suppressed at higher fields. It appears that the ratthng 
frequencies are too high, in comparison with the ESR fre- 
quencies, and the spatial inhomogeneities are completely 
averaged out. 

We ascribe the lack of clear signatures of the rattling 
behavior in the Q-band ESR spectra to the relatively high 
fields {H « 1 T) used in the experiment. As observed 
in ESR experiments in a true spin-glass systenpSl^ higher 
fields tend to quench the random anisotropy contribu- 
tion. Our claim is that a field as high as IT is enough 
to suppress the small anisotropics implied by the Eu^+ 
dynamical behavior thus suppressing the broadening of 
the spectra when Tq is approached. 

Figure |4] shows that no comparable phenomenon takes 
place in the clathrate compound. Figure [IF a) shows that 
at low-T, the spectra are anisotropic due to the prox- 
imity of the magnetic transition (Tc = 35 K^^l) . Above 
this low-T region (T > 70 K), the spectra are isotropic 
with g = 1.999(3) and broaden in a Korringa like relax- 
ation. Figure |4Fb) evidences the relatively narrow ESR 
line when compared to the skutterudite compounds. As 
the transition T is approached, and below this temper- 
ature, AH broadens rapidly and continuously. The g- 
values also refiect clearly the low-T anisotropy and mag- 
netic transition. 

No clear connection with the Eu^'^ rattling modes is 
observed in this experiment. Given the specific features 
of the behavior of the filler atom in this clathrate {off- 
center rattling), we would expect more signatures aris- 
ing from crystal field inhomogeneities in this experiment 
than in the experiments with the skutterudites. This, 
and the Q-band measurements, lead us to rule out crys- 
tal field effects as the origin of the phenomena observed 
in figures flla)-(b) and[2Fa)-(b). Both results considered 
together, give strength for our claim that the unique be- 
havior of the ESR spectra of the skutterudites are due to 
the coupling between the rattling modes and electronic 
degrees of freedom. 

Although we observe a metallic relaxation, the pro- 
cess is slow {b = 0.86 Oe/K), which agrees with previ- 
ous results pointing that EugGaigGeao is a poor metal- 
lic systeniiSl This Korringa rate is compatible with the 
measured Ag = 0.006(5), pointing to a q-independent 
Jfs constant and absence of multi-band effects. 

It is noteworthy that AH in the skutterudites and in 
the clathrates are of the same order of magnitude at low- 
T (see figures [4] and [3] around T w 50 K) . In general, 
AH in clean concentrated metallic systems will be de- 
termined by strong dipolar interactions, which broaden 
the line, and spin-spin Exchange narrowing effects. Since 
the Eu^"*" ions are slightly more apart from each other in 
the clathrates, one would expect dipolar interaction to be 
weaker in the clathrates and hence AH would be smaller, 
as it is verified. Since the latter two effects are temper- 
ature independent, the significant difference in AH at 



high-T should be ascribed to the different Korringa rates 
of the Korringa-like relaxation^^, discussed above. The 
exact nature of a relaxation process in a concentrated 
metallic system is very hard to determine. However, it 
should be related to an Exchange interaction with con- 
duction electrons, and that is why we are carefuU to refer 
to the linear broadening of AH as a "Korringa-like" pro- 
cess. All the reported material parameters relevant to 
this coupling, su ch as the density of electronic states at 
the Fermi surface^SESl^ favors the idea that this coupling 
(and therefore the relaxation) should be stronger in the 
skutterudites. 

Contrary to the skutterudites studied here, we note 
that for this compound, Raman scattering studiea^S have 
shown that the Eu^'^ rattling energy {0e ~ 25 K) is 
lower than that of the magnetic ordering {Tc = 36 K). 
Hence, a spin-glass like behavior below 6e may also be 
prevented from occurring by the magnetically ordered 
state. Furthermore, the itinerant d orbital contribution 
to the Fermi surface, in the case of the skutterudites, 
is most likely originated from the Ru (Os) d atomic or- 
bitals, indicating that the skutterudite cage is somehow 
stiffer (more correlated) than the clathrate cage. Indeed, 
the Ag analysis point for a greater coupling of localized 
and itinerant states in the skutterudites, which should 
be important to correlate the rattling behavior with the 
Eu-Eu interaction. 

In some skutterudites, there are earlier reports for tem- 
perature dependence of the valence state of the guest i on, 
as was reported for the related compound YbFe4Sbij2ll. 
This result raises the question about the valence state 
of the Eu ion on the compounds here investigated. We 
should point out, however, that the issue of the Yb va- 
lence state in YbFe4Sbi2 was subsequently revised in the 
literature's determining that the Yb ions are in a stable 
divalent state. In our compounds, no significant frac- 
tion of Eu^+ was found for EuT4Sbi2 (T =Ru, Os) nor 
it was found any sign for temperature dependency of 
the Eu^+ fraction in these materials^. In EuFe4Sbi2, 
a fraction of 10 — 15 % of Eu"^"*" was reported. How- 
ever, again, no sign for temperature dependency of this 
fraction was found, e ven in detailed studies of X-ray ab- 
sorption spectroscopjUSHZl, xhe insets of figures k[ a)- (b) 
show a Curie-like behavior for the ESR intensity, that 
should also be taken as another piece of evidence for a 
nearly stable Eu^+ configuration. 



IV. CONCLUSION 

In conclusion, we have provided evidence, from Eu^^" 
X-band ESR measurements, for inhomogeneous spin dy- 
namics of the Eu^+ spins in the EuT4Sbi2 (T= Ru, Os) 
skutterudites, triggered by the Eu^"*" rattling modes. Our 
Q-band measurements were discussed in terms of present- 
ing evidence for q-dependence of J{q) exchange coupling 
and for covalent mixing between localized / and itinerant 
d orbital states. These two findings were related, respec- 



tively, to the inhomogeneous spin dynamics of the Eu^"*" 
spins and as an evidence for a stronger coupHng between 
localized and itinerant states in the skutterudites than in 
the clathrate. 

Our picture is that the behavior of AH in the X-band 
ESR measurements is due to an electron-phonon coupling 
between the rattling phonon modes and the electronic de- 
grees of freedom. In regions I and II of figures p[a) and 
[2[b), the itinerant d orbital electronic states are coupled 
with oscillating and non-oscillating 4/ electronic states, 
which gives rise to random anisotropics of the exchange 
coupling and subsequently to the inhomogeneous spin dy- 
namics found in these regions. 

In analogy with the phonon-glass scenario discussed for 
cage systems, we discussed our findings in the framework 
of a spin-glass type of spin dynamics. In this sense, high 
fields tend to suppress the ESR line broadening when Tq 
is approached, by quenching the random anisotropics of 



these systems. We suggested that a similar effect takes 
place in our experiment, preventing the observation of 
inhomogeneous dynamics of the Eu^+ spins in Q-band. 
Further theoretical and experimental investigation may 
unravel the exact characteristics of this spin-glass like 
state and also explore the perspective of other phenom- 
ena emerging from the interaction between conduction 
electrons and the rattling phonon modes. 
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